Complete QCD Corrections to Hadronic Z-Decays in Order a A s 
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Corrections of order ai for the axial singlet contributions for the decay rate of the Z-boson into 
hadrons are evaluated in the limit of the heavy top quark mass. Combined with recently finished 
(D(ai ) calculations of the non-singlet corrections, the new results directly lead us to the first complete 
0(a*) prediction for the total hadronic decay rate of the Z-boson. 

The new 0(o4) term in Z-decay rate lead to a significant stabilization of the perturbative series, 
to a reduction of the theory uncertainty in the strong coupling constant a 3 , as extracted from these 
measurements, and to a small shift of the central value. 
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The precise determination of the Z-boson decay rate 
into hadrons at LEP [l[ has led to one of the most precise 
determinations of the strong coupling constant a s . From 
the experimental side, in view of the fully inclusive nature 
of this measurement, the result is fairly robust, in partic- 
ular since it is insensitive to simulations of the hadronic 
final state. Hence the error is essentially dominated by 
the statistical uncertainty. From the theory side the ad- 
vantage of the measurement is its high energy, and as a 
result, the irrelevance of nonperturbative and power-law 
suppressed terms. The smallncss of a s at high energies 
then leads to a rapid decrease of higher order corrections 
in the perturbative series and, correspondingly, to a sig- 
nificant reduction of the theory error. 

A variety of methods has been suggested to estimate 
the remaining uncertainty in the theory prediction. Us- 
ing the last calculated term is probably the most con- 
servative approach, varying the renormalization scale /i 
within an energy range characteristic for the problem 
(e.g. Mz/3 < n <3Mz) is frequently used, albeit with 
considerable ambiguity in the actual choice of the region 
of the /x-variation. In order to reduce the theoretical un- 
certainty in the extraction of a s to a level significantly 
smaller than the experimental one (which amounts to 
±0.0026 at present [l|]), the knowledge of the corrections 
of 0(ctg) is necessary. At the same time this calcula- 
tion opens the window for a considerable improvement 
in the a s -determination at Giga Z , the project of a high- 
luminosity linear collider operating at the Z-resonance 
(see e.g. Q, where a precision of 0.0005 to 0.0007 
has been advertised). The dominant part of the a*- 
corrections, the "non-singlet" -piece, has been evaluated 
in 0] . This has lead to a slight shift of the central value 
oia s upwards from 0.1 185 ±0.0026 to 0.1190±0.0026 Q 
and a reduction of the theory error far below the error 
of 0.0026 from experiment. However, as noted already in 
Jjjj, for a complete evaluation of the decay rate in 0(a 4 s ) 
an additional set of corrections, namely those for the 



"singlet" contributions, is required. For the axial cur- 
rent correlator these start at 0(a 2 ) [3, [f|, for the vector 
correlator at O(a^). Both of them are presently known 
to third order in a s only Hence, for a completely 

consistent C(a*) extraction of the strong coupling the 
extension of these results by one order in a s is required. 

Before describing this calculation in detail, let us 
briefly recall the basic structure of QCD corrections to 
the correlator of the electromagnetic and the neutral cur- 
rent, respectively, their similarities and their main differ- 
ences. After splitting off inessential kinematic factors, 
the absorptive part of the current-current correlator of 
the electromagnetic current is expressed by the familiar 
i?-ratio 



R cm = 3 



E^ns + (E<?/) 2 ^ 

/ / 



(1) 



where rj^ s and r$ stand for the (numerically dominant) 
non-singlet and the singlet part respectively. The corre- 
sponding decomposition for the correlator of the neutral 
current involves the following four terms 
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2If and sw defined as 
Here all but the top quark 



with Vf = 2If — iqfSyyi o>f 
effective weak mixing angle, 
are assumed to be massless. 

(Mass corrections to both vector and axial vector cor- 
relator due to other massive quarks are dominated by the 
bottom quark and can be classified by orders in m 2 /M z 
and a s . Up to 0(c^ m 2/M§) and 0{a 1 s m\IM\) they can 
be found in [ll[, as well terms of order o? s m\\jM\ (const 
+ log ml/Ml) and a 2 s ml/M? (const + log mg/M|) 
that arise from axial vector singlet contributions. Terms 
of order o? a m\jM z and afm^/Mf can be found in [lit ] 
and [IH respectively. Corrections of order a 2 s m 2 z /m 2 and 
a l m z/ Tn t f rom singlet and non-singlet terms are known 
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from 0, [H, [3] and [l5| respectively. These are important 
for the actual a s -determination, but will not be discussed 
further in the present paper.) 

From the prefactors of the non-singlet contributions 
in electromagnetic, vector and axial correlator it is ev- 
ident that different quark flavours contribute incoher- 
ently, hence additive to the rate. Thus their contribution 
is significantly enhanced in comparison with the singlet 
terms where amplitudes from different flavours interfere 
destructively, with prefactors (J2j qj) 2 and (Y^,f v f) 2 f° r 
the electromagnetic and neutral current respectively. 

Non-singlet contributions are present at the parton 
level and the QCD corrections are known in second flij . 
third 6, 7] and fourth Q order in a s . In terms of Feyn- 
man diagrams, non-singlet contributions are character- 
ized by the fact that one quark loop connects the two ex- 
ternal currents (Fig. [Taj) . In the absorptive part of this 
fcrmion loop no top quark is present due to kinematic 
reasons, whence the non-singlet functions are identical 

r NS = r NS = r NS- 

In the case of singlet contributions of the vector current 
the two currents couple to two different quark loops (Fig. 
Tb|) requiring a three-gluon intermediate state. Corre- 
spondingly the leading term is of 0(a 3 ) and has been 
obtained already long time ago [1, Q ■ The NLO correc- 
tions to this result are of 0(a A s ). They serve to soften the 
strong scale dependence of the 0(a 3 ) result, stabilize the 
theory prediction and will be the subject of this paper. 

The situation is different in the case of the singlet ax- 
ial vector current correlator. The axial couplings of the 
two members of an isospin doublet are opposite equal. 
Hence their singlet contributions vanishes, if the corre- 
sponding quark masses are equal. This approximation 
is valid for the two lightest quark doublets. The only 
remaining contribution originates from the combination 
of bottom and top quarks with their specific mass hier- 
archy m 2 « M§ << m 2 (Fig. [Tc]) . In this case the 
contribution starts at 0(a 2 ) and is further enhanced by 
the "large" logarithm logjirif/M 2 ,) [J , la- Corrections of 



0{a 3 s ) have been calculated in |8l4lfj|. those of O(a^) will 
be the subject of this paper. 

The evaluation of the NLO terms of r v requires the 
calculation of the absorptive parts of five-loop diagrams 
with massless propagators which, with the help of some 
complicated combinatorics based on the -Reoperation 
[l7j |. can be boiled down to the calculation of four- loop 




FIG. 1: Different contributions to r-ratios: (a) non-singlet, 
(b) vector singlet and (c) axial vector singlet. 



propagator diagrams. The latter have been computed 
via reduction to 28 master integrals, based on evaluating 
sufficiently many terms of the 1 / D expansion flij of the 
corresponding coefficient functions [I9j |. This direct pro- 
cedure required huge computing resources and was per- 
formed using a parallel version [2fJ of FORM 2J| . The 
master integrals are reliably known from [22T - [2~ij . The de- 
tails of the calculation, the results in analytic form and 
their relation to the Gross-Llewellyn Smith sum rule will 
be given in [27j |. 

The evaluation of the NNLO terms of i?g. t b involves 
again absorptive parts of five-loop diagrams with mass- 
less propagators, however, in addition also absorptive 
parts of four-loop diagrams combined with one-loop mas- 
sive tadpoles, etc. down to one-loop massless diagrams 
together with four-loop massive tadpoles. The latter have 
been computed with the help of the Laporta algorithm 
25] implemented in Crusher [26[ . The methods employed 
in our calculations, together with the results will be de- 
scribed in more detail in 



The result is valid in the limit 



M 2 



< 4 Mr, an ex- 



cellent approximation as evident from the lower orders. 
The relative importance of the various terms is best seen 
from the results for the various r-ratios introduced above, 
expressed in numerical form 



r NS =1 + a s + 1.4092 a 2 s - 12.7671 a 3 s - 79.9806 a* 
rX = - 0.4132 a 3 - 4.9841 at , 



r&,b= (-3-0833 + l t )a 2 s 



(3) 



+ (-15.9877 + 3.7222 l t + 1. 9167 if) a 3 s 
(49.0309 - 17.6637 l t + 14.6597 1 2 + 3.6736 I 3 ) a 
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with a s = a s (M z )/ir and l t = ln(M|/M t 2 ). Since all 
three r-ratios are separately scale invariant, the corre- 
sponding results for a generic value of a s (fi) can easily 
be reconstructed. Using for the pole mass M t the value 
172 GeV the axial singlet contribution is given in numer- 
ical form by 



r S;t,b 



= -4.3524 a 2 - 17.6245 a 3 + 87.5520 at 



(4) 



Collecting now all QCD terms, the decay rate of the Z- 
boson into hadrons can be cast into the following form 

r z = r o i?- = |0|i? nc . (5) 

Here all electroweak corrections are assumed to be col- 
lected in the prefactor Tq, and the forementioned mass 
corrections are ignored as well as electroweak and mixed 
QCD-electroweak corrections 23-30]. Thus the i?-ratio 
is now known up to 0(a^) 

R nc =20.1945 + 20.1945 a s 

+ (28.4587- 13.0575 + 0) a 2 

+ (-257.825 - 52.8736 - 2.12068) a 3 

+ (-1615.17 + 262.656 - 25.5814) a 4 s , (6) 
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with 



Mr 



0.231. The three terms in the brackets dis- 



play separately non-singlet, axial singlet and vector sin- 
glet contributions. 

Let us now evaluate the impact of the newly calcu- 
lated terms on the a s -determination from Z-decays. Fol- 
lowing our approach for the non-singlet terms (where 
a shift Sa s — 0.0005 had been obtained consis- 
tent with an analysis }3l| based on results of the elec- 
troweak working grou p | l|] and a modified interface to 
ZFITTER v. 6.42 [II 111] and confirmed by the G-fitter 
collaboration [32,30,31]), we consider the quantity i? nc 
as "pseudo-observable". With a starting value R nc — 
20.9612, if evaluated for a s = 0.1190 and without the a\ 
singlet terms, a shift 8a s = —0.00008 is obtained after 
including the newly calculated contributions. 

As discussed in [3J, the non-singlet a* term leads to a 
considerable stabilization of the theory prediction, and, 
correspondingly, to a reduction of the theory error. A 
similar statement holds true for the singlet contribution. 
To illustrate this aspect, the dependence on the renor- 
malization scale fj, is shown in Fig. [5] for tns , tY and 
7g. t b . The relative variation is significantly reduced in 
all three cases. In particular for the vector singlet case 
we observe a shift of the result by about a factor 1.45 
(for fi — Mz) and a considerable flattening of the result. 
Using for example the Principle of Minimal Sensitivity 
(PMS) [35| as a guidance for the proper choice of scale, 
/i = 0.3 Mz seems to be favoured, leading to an amplifi- 
cation of the LO result by a factor 1.68 (if the latter is 
evaluated for /i = Mz, as done traditionally). 

Let us assume that the remaining theory uncertainties 



from tns , rY and . 



S;t,b 



can be estimated by varying fj, be- 



tween M^/3 and 3Mz and using the maximal variation 
as twice the uncertainty Sr. This leads to ST^s = 0.101 
MeV, STY. = 0.0027 MeV and STg = 0.042 MeV. Even 
adding these terms linearly, they are far below the exper- 
imental error of ST exp = 2.0 MeV [36| , In combination 
with the quadratic and quartic mass terms, which are 
known to 0(a\) and 0(af) respectively, this analysis 
completes the QCD corrections to the Z decay rate. 

Let us also comment on the impact of the a s singlet 
result on the measurement of i? cm at low energies, i.e. in 
the region accessible at BESS or at B-factories, say be- 
tween 3 GeV and 10 GeV. Considering the large luminosi- 
ties collected at these machines, a precise a s determina- 
tion from i? em seems possible [38( . In the low energy re- 
gion only rg and r^ s contribute. Since X)/= u d s If = 0j 
the singlet contribution vanishes in the three flavour case. 
If we consider the region above charm and below bottom 
threshold, say at 10 GeV, only u,d,s and c quarks con- 
tribute, the relative weight of the rg in eq. (p} is given by 
(H Qf) 2 /(J2 Qf) — 2/5, and thus is fairly suppressed. At 
energy of 10 GeV, in the absence of open bottom quark 
contribution, it seems appropriate to analyze the results 



in an effective four flavour theory with 
rY. = -0.41318 a 3 s {n) - (5.1757 + 2. 58241n^ 2 /s) 

As shown in Fig. [3l it is evident that the scale depen- 
dence is softened in NLO. Again a scale /i around 0.3 y/s 
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FIG. 2: Scale dependence of (a) non-singlet r^s, (b) vector 
singlet rg and (c) axial vector singlet rg. t b . Dotted, dash- 
dotted, dashed and solid curves refer to 0(a a ) up to OipQ) 
predictions. a s (Mz) = 0.1190 and m = 5 is adopted in all 
these curves. 
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minimizes the NLO corrections. 




FIG. 3: Scale dependence of the vector singlet r s around 
10 GeV. Dashed and solid curves refer to 0{a 3 s ) and 0(a$) 
predictions, ni — 4 and a 3 (10GeV) = 0.1806 as obtained 
with the use of package RunDec [37| have been assumed. 

In conclusion we want to mention that all our calcu- 
lations have been performed on a SGI ALTIX 24-node 
IB-interconnected cluster of 8-cores Xeon computers us- 
ing parallel MPI-based as well as thread-based [39} 
versions of FORM [2l[ . For evaluation of color factors we 
have used the FORM program COLOR Ep. The dia- 
grams have been generated with QGRAF [41[. This work 
was supported by the Deutsche Forschungsgcmeinschaft 
in the Sonderforschungsbereich/Transregio SFB/TR-9 
"Computational Particle Physics" , by Graduiertenkolleg 
1694 "Elementarteilchenphysik bei hbchster Energie und 
hochster Prazision" and by RFBR grants 11-02-01196 
and 10-02-00525. 
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